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ABSTRACT 

The  refractive  index  structure  coefficient,  CR,  is  an  important 
parameter  for  optical  propagation  in  the  atmosphere.  A  theoretical  dis¬ 
cussion  of  structure  functions  is  presented  to  show  the  origin  of  the 
refractive  index  and  temperature  structure  functions.  The  equipment  used 
to  measure  temperature  fluctuations  in  the  atmosphere  is  discussed,  and 
the  calculation  of  the  refractive  index  structure  coefficient  from  these 
thermal  measurements  is  shown.  The  details  of  the  experimental  evaluation 

of  C  and  the  results  of  C  calculations  are  also  included, 
n  n 
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INTRODUCTION 


In  a  turbulent  medium  the  velocity  of  the  fluid  flow  et  any  point  ie 
defined  by  a  vector  which  varies  with  the  turbulence.  The  turbulence  can 
be  characterized  by  nine  structure  functions  relating  the  velocity  vectors 
at  tvo  points  in  the  flow.  A  measure  of  the  homogeneity  of  the  medium  can 
be  expressed  by  the  structure  coefficients  derived  from  these  functions. 

Using  a  theoretical  model  developed  by  V.  I.  Tatarski  for  atmospheric 
turbulence  affecting  wptical  propagation,  a  method  of  directly  measuring 
the  temperature  fluctuations  between  two  points  normal  to  the  wind  veloc¬ 
ity  is  employed  to  evaluate  the  refractive  index  structure  coefficient  0  . 
When  this  meteorological  determination  of  CR  is  correlated  with  the 
observed  characteristics  of  monochromatic  light  propagated  through  the 
atmosphere,  the  applicability  of  the  body  of  theory  developed  by  Tatarski 
and  others  to  optical  propagation  in  the  atmosphere  can  be  determined. 

THEORETICAL  BACKGROUND 

The  structure  functions  of  a  turbulent  flow  are  of  the  form 


D^flr)  =  v1(r1)  -  V^}]  '  VV]»  f°r  i,J  =  1'2'3'  (1'1) 


where  the  V  's  are  the  velocity  components  of  the  medium  at  the  points 

-A  -A  Jk  -A 

and  and  r  =  -  r^.  The  bar  over  the  right-hand  side  of  Equation 

(l.l)  indicates  time  averaging  of  the  product  beneath  it.  (The  time 
average  gives  a  measure  of  the  expected  value  of  the  quantity.)  Thus,  the 
structure  functions  are  related  to  the  expected  values  of  the  covariance 
of  the  velocities  at  the  two  points. 

If  a  condition  of  local  isotropy  is  assumed  to  exist  in  the  turbulent 
medium  and  the  two  points  and  are  in  a  plane  perpendicular  to  the 
flow  end  within  the  range  of  iso  .ropy,  the  tensor  D. .  can  he  expressed  as 
a  function  of  r  =  |r|j  i.e.,  the  distance  between  the  two  points.  The 
longitudinal  structure  function  I-r(r)  and  the  transverse  structure  func¬ 
tion  D^(r)  then  determine  ^(r).  Making  the  further  assumption  of 
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homogeneity  of  the  flov,  one  can  express  these  two  structure  functions  in 
tern*  of  each  ether  and  thus  characterize  the  turbulent  flow  by  one  struc¬ 


ture  function  depending  only  on  r.  When  the  separation  of  the  two  points 
ia  large  compered  to  the  inner  scale  of  the  turbulence,  denoted  by  1 

t  o 

(the  dimension e  of  the  smallest  turbulent  eddies  in  the  flov)  but  small 

compered  with  the  outer  scale,  Lq  (the  largest  dimension  for  which  the 

medium  can  be  assumed  to  be  isotropic),  tbe  structure  functions  are 

directly  proportional  to  r2/5.  The  validity  of  this  "two-thirds  law," 

14  A 

first  obtained  by  Kolmogorov  and  CJbukov,  has  been  investigated  by  other 
experimenters  ’  and  found  to  hold  reasonably  well.  The  square  root  of 
the  coefficient  of  proportionality  is  called  the  structure  coefficient  of 
the  flow,  and  the  size  of  this  coefficient  is  a  measure  of  the  local 


inhomogeneity  of  the  medium.  For  a  steedy-state  laminar  flow,  this  eoef- 
ficlmt  would  be  zero  as  the  medium  would  be  totally  homogeneous  and  the 
velocities  a*  all  points  in  the  medium  the  Bame. 


Treating  temperature  as  a  "conservative  passive  Bdditive"  of  the 
atmosphere,  l.e.,  a  parameter  that  has  &  direct  influence  on  atmospheric 
turbulence,  Obuhov^  related  the  structure  functions  of  turbulent  flow  to 
the  structure  function  of  the  atmospheric  temperature  field.  From  this 
work  Tatars Ki  derived  an  expression  for  the  refractive  Index  structure 
function: 


Dn(rJ  "  Cnr 2/3  '  for  !©  «  r  «  I-o  >  (1-2) 

where  CQ  is  tbe  structure  coefficient  of  atmospheric  refractive  index  and 
may  be  thought  of  as  a  measure  of  the  optical  "strength"  of  atmospheric 
turbulence. 

The  analogous  expression  for  the  temperature  structure  function 
developed  by  Cfeukov  is 

DT(r)  -  C T2r2/5  ,  for  1Q  «  r  «  Lq  ,  (1.3) 

1 - 

Superscript  nmbere  denote  referenoea  vkieh  may  be  found  on  page  IB. 
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where  is  the  structure  coefficient  of  temperature.  Equations  (1.2) 
end  (1.3)  give  the  theoretical  relationships  necessary  to  evaluate  Cn  from 
the  measurement  of  atmospheric  temperature  fluctuations. 


INSTRUMENTATION 

The  Instrument  used  to  measure  the  temperature  fluctuations  between 
two  points  in  the  atmosphere  is  essentially  s  Flow  Corporation  Hot-Wire 
Anemometer  modified  for  use  ss  a  high-speed  atmospheric  thermometer.  The 
Flow  Corporation  equipment  consists  of  two  high-speed  temperature  probes, 
the  Anemometer  Bridge  Model  HWB-3,  a  secondary  voltage  supply  Model  HWI-3, 
and  a  Random  Signal  Voltmeter  Model  L2A-1,  Figure  1.  Data  are  read  with 
a  digital  voltmeter  and  recorded  on  paper  tape  by  a  digital  clock/printer. 
The  digital  voltmeter  samples  the  output  signal  of  the  continuous  reading 
Random  Signal  Voltmeter  at  1-second  to  1-minute  intervals.  B»e  clock/ 
printer  records  the  time  and  the  voltage  reading  for  each  sample. 

The  temperature  probes  utilize  fine  tungsten  wires  of  10  microns  in 
diameter  and  approximately  0.3  centimeter  in  length  as  sensing  elements. 
For  precision  temperature  measurements,  an  operating  current  of  1  ma  is 
used  in  these  vires.  The  resistance  of  the  tungsten  wires  In  the  range 
of  atmospheric  temperatures  is  essentlnlly  a  function  linear  with  tempera¬ 
ture  and  may  be  written 

T  *  KpR  +  TR  .  (2.1) 


T  is  in  degrees  Kelvin,  and  R  is  in  ohms.  TR  is  a  constant  equal  to  the 
temperature  corresponding  to  zero  resistance  on  the  straight  line  curve. 
When  taking  the  temperature  difference  measurement  between  the  two  parches, 
a  difference  in  the  resistance  constant  for  different  probes  Introduces 
an  error  in  the  resistance  measurement.  The  probes  are  balanced  to  give 
a  zero  temperature  difference  for  any  change  in  average  air  temperature 
in  order  to  minimise  this  error.  The  difference  in  the  two  resistance 
constants  is  small,  and  therefore  the  error  made  in  measuring  the  small 
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&  c* 


short-time  temperature  differences  relating  to  CQ  Is  negligible.  The 
average  of  the  tvo  constants  is  used  for  calculation  proposes.  The  dif¬ 
ference  in  temperature  between  the  tvo  probes  can  be  written  approximately 
as  the  difference  in  their  resistance  multiplied  by  the  average  resistance 
constant  K^, 

AT  -  Kpflti  .  (£.2) 

Vlth  a  current  of  1  ms  in  the  vires,  resistance  in  ohms  is  directly  equiv¬ 
alent  to  voltage  measured  In  millivolts.  Therefore, 

AT  -  KpAV  .  (2.3) 

The  thin  tungsten  vires  Insure  fast  response  of  the  probes  to  atmo¬ 
spheric  temperature  fluctuations,  but  the  large  surface  area  to  volume 
ratio  makes  the  probes  sensitive  to  solar  heating.  It  is  assumed  any 
solar  heating  is  equal  for  both  probes  and  so  will  have  no  effect  on  the 
temperature  difference  measurement. 

The  two  probe  voltages  are  balanced  through  the  bridge  circuits  of 
the  HWB-3  equipment  using  the  galvanometer  Included  In  this  unit  until  a 
null  la  achieved  across  the  bridge.  This  calibration  is  made  with  the 
temperature  probes  enclosed  in  a  reflecting  wooden  box  to  insure  that 
both  probes  ere  exposed  to  the  same  average  air  temperature. 

After  the  tvo  probe  voltages  are  nulled,  the  probes  are  securely 
mounted  at  a  separation  distance  of  30  centimeters  horizontally  about  the 
axis  of  a  wind  vane.  Figure  2.  The  wind  vane  maintains  the  probes  In  a 
perpendicular  orientation  with  respect  to  the  wind  flow.  The  30-centimeter 
distance  satisfies  the  conditions  of  Equations  (1.2)  and  (1.3)  as  the 
inner  scale  of  atmospheric  turbulence  is  on  the  order  of  a  few  millimeters, 
while  the  outer  scale  affecting  optical  propagation  Is  on  the  order  of  2 
meters,  the  beam  height  above  the  ground.  When  the  two  probes  are  exposed 
to  the  air,  the  fluctuations  in  atmospheric  temperature  produce  a  continuous 
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Tenperature  probes  mounted  on  wind  vane 


voltage  signal  in  the  bridge  circuit  corresponding  to  the  differences  in 
temperature  at  the  two  probes  as  indicated  in  Equation  (2.3).  This  signal 
is  then  amplified  and  fed  into  the  12A-1  Random  Signal  Voltmeter. 

The  Random  Signal  Voltmeter  processes  voltages  by  means  of  an  inter¬ 
nal  thermal  block  with  a  time  constant  of  approximately  16  seconds.  The 
response  of  the  thermal  block  varies  as  the  square  of  the  input  voltages. 
Thus,  the  Random  Signal  Voltmeter  reading  is  the  16-second  time  average 
of  the  square  of  the  voltage  differences  and  1b  directly  proportional  to 
the  time  averaged  square  of  the  temperature  gradient  at  the  two  probes  by 

-  o 

the  constant  ; 

dtf2  «=  K^AV2  .  (2. h) 

Taking  the  square  root  of  Equation  (2.k)  to  obtain  the  time  average  root 
mean  square  temperature  gives 

“ms  -  y*ms  ■  (2-5> 

A  potential  source  of  error  arises  from  the  small  signal  to  noise 
ratios  encountered  in  the  measurement.  The  use  of  internal  battery  volt¬ 
age  supplies  for  probe  biasing  eliminates  the  possibility  of  introducing 
interior  random  or  cyclic  noise.  The  remaining  noise  is  due  to  random 
signal  pickup  in  the  cables  leading  from  the  probes  to  the  bridge.  As 
this  noise  level  (read  during  calibration)  is  relatively  constant,  it  can 
be  subtracted  from  the  final  signal  to  greatly  reduce  this  error.  Appen¬ 
dix  A  gives  the  detailed  calculations  of  the  temperature  difference  includ¬ 
ing  the  correction  for  noise  level. 

DETERMINATION  OF  C 

n 

From  the  preceding  discussion  of  the  structure  functions  of  a  turbu¬ 
lent  medium,  ve  can  define  the  structure  functions  of  temperature  and 
refractive  index  as 
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(3.1) 


2 

Dj^r)  *  JoKi^)  -  T(r2)J 

Dn(r)  =  Jn^)  -  n(r2)j  (3.2) 

where  r  -  ji^  -  r2| 

and  the  average  bars  indicate  time  averaging.  Equation  (3.1)  is  an 
expression  for  the  time  averaged  Bquare  of  the  temperature  difference 
between  r^  and  r£  and  can  be  related  to  the  readings  of  the  Random  Signal 
Voltmeter  by  Equation  (2.4) 

DT(r)  =  dff2  =  K^AV2  .  (3.3) 

Using  Equations  (1.2)  and.  (1.3)  presented  earlier,  we  obtain 

Cn  "  CT2(Dn/DT)  *  (5A) 

The  quantity  D  /D_  in  this  expression  is  a  function  of  temperature,  but 
n  T  ii 

the  dependence  is  as  l/T  ,  where  T  is  in  degrees  Kelvin.  For  small  tem¬ 
perature  changes  D^/D^  can  be  considered  a  constant.  Appendix  B  justifies 
this  approximation  and  gives  the  method  used  to  evaluate 

Using  Equations  (1.3)  and  (3-3),  Equation  (3-4)  can  be  written  as 

CQ2  =  £s^t'2^  •  (5*5) 

Taking  square  roots  and  combining  the  two  constants  gives  a  concise 
expression  for  the  structure  coefficient 

CQ  =  CV^gr"1/3  ,  (3.6) 

where  the  constant  C  is  equal  to  K^D^/D^)  '  . 
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The  refractive  index  structure  coefficient  C  was  evaluated  for 

n 

18  hours  on  18  August  1967  from  0600  to  2200  hours.  The  results  of  the 
measurements  are  plotted  against  time  of  day  in  Figure  3.  Rie  Random 
Signal  Voltmeter  readings  were  averaged  to  10-minute  intervals  for  this 
analysis.  The  isolated  points  for  the  first  few  readings  result  from  gaps 
in  the  data  caused  by  repeated  calibration  checks  at  the  start  of  the 
experiment.  Hie  maximum  value  of  is  indicated  at  1300  hours,  ESOT, 
and  therefore  coincides  with  the  sun’s  zenith. 

Data  were  taken  on  22  August  1967  from  0900  hours  to  iVOO  hours. 

The  results  of  the  calculations  for  this  day  are  plotted  in  Figure  k 
(5 -minute  average)  and  Figure  5  (lO-minute  average).  Hie  maximum  value 

-ft 

for  Cn  again  occurs  near  1300  hours.  Also  important  is  the  full  range 
of  Cn  on  both  these  days  exhibiting  variations  of  nearly  100  to  1. 

It  is  paradoxical  that  a  structure  function  describing  an  optical 
parameter  is  better  quantified  by  a  thermal  rather  than  an  optical  mea¬ 
surement.  However,  it  has  been  demonstrated  that  the  index  structure 

coefficient,  C  ,  cannot  be  correctly  computed  through  the  well-known 
n  ** 

equation  given  by  Tatarski  for  the  amplitude  variance  in  terms  of  the 

7 

path  length,  wave  number,  and  the  index  structure  coefficient.  It  is 
apparent  that  measurements  of  appropriate  meteorological  parameters  are 
needed  before  an  adequate  understanding  of  optical  propagation  in  the 
atmosphere  can  be  gained. 


It  is  expected  that  Cn  should  reach  a  peak  at  this  time  as  the  solar 
flux  to  the  earth's  surface  is  at  a  maximm  producing  the  greatest  heat 
transfer  at  the  air-ground  interface. 

«* 

See  Reference  0,  p.  188 ,  Equation  (9.43). 
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TIME  AVERAGED  BY 


minute  average 


22  AUGUST  1967  (5-MINUTE  A  VERA  GO 


HOURS 

versus  hou  of  day  (EDST),  22  August  1967,  5-minute  average 


22  AUGUST  1967  (10-MINUTE  AVERAGE) 


HOURS 

versus  hour  of  day  (EfST),  ZZ  August  1967,  10-minute  average 
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APPENDIX  A 


The  voltage  at  the  meter  output  of  the  Flow  Corporation  12A1  is 
related  to  the  reading  on  the  RMS  meter  in  the  following  way: 


E 


mo 


(A-l) 


where  E  -  voltage  at  meter  output  terminals 
no 

Eln  ■  meter  reading  (RMS  value  of  input  voltage) 

Efg  *  full-scale  voltage  on  the  range  in  use 
K  *  a  meter  constant,  approximately  5  mV. 

The  thermal  probee  have  the  following  linear  relationship  with 
temperature: 

T  =  KjjR  +  T^  ,  ( A-2  )• 

where  T  is  the  temperature  (degrees  Kelvin),  R  is  the  resistance  (ohms), 
and  Kp  and 

The  temperature  difference  then  is 

AT  =  KpflR  .  (A- 5) 


TL  are  constants  that  are  found  by  calibrating  the  probes. 


If  the  probes  are  biased  with  a  constant  current  of  1  ma, 

AT  *  K^AV  ,  (A-4) 

where  AV  is  measured  in  millivolts. 

If  both  probes  were  to  have  the  same  Kp,  then 

%hs  ’  Wne  ■ 

where  AV  is  the  voltage  (divided  by  the  appropriate  gain  factor)  displayed 
on  the  random  signal  voltmeter.  Usually  the  two  K^'s  are  nearly  the  sane; 
therefore,  averaging  the  slopes  of  the  two  probe  curves,  the  approximation 

Ip  «  (A-6) 
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&■  UNd. 


Relating  Bquatlons  ( \-l)  and  (A-5)  we  now  find 


(A-7) 


The  factor  of  1/50  la  included  because  we  are  aasuaing  all  meaeuremeuta 
will  be  aede  with  galvanoawter  range  switch  on  the  0.2  mV  scale. 

Situation  (A-7)  would  be  correct  If  the  noise  level  of  the  Instrument 
were  aero.  A  noise  level  is  present,  however,  and  oust  be  subtracted 
fraa  the  readings  to  obtain  an  accurate  temperature  difference  aeasure- 
aent.  Accounting  for  the  system  noise  gives  finally, 


Off, 


IMS 


'  S«o®fs. 


\l/2 


J 


(A-8) 


where  the  *n”  subscript  denotes  noise  readings  recorded  during  calibration 
of  the  equipment. 


AFPEHDIX  B 


In  order  to  show  that  the  quantity  Dq/Bt  is  proportional  to  1/T  , 
where  T  is  temperature  in  degrees  Kelvin,  the  following  relationship  can 
be  used  giving  the  refractive  index  of  air  as  a  function  of  temperature, 
pressure,  and  wavelength: 

(n  -  1)  -  77.6P/T( 1  +  0.0075X2)  x  10'6  (B-l) 

fop  pressure  P  in  millibars  and  wavelength  X  in  meters.  For  small  tem¬ 
perature  changes  atmospheric  pressure  can  be  taken  to  be  constant  allowing 
the  change  in  index  of  refraction  with  respect  to  temperature  to  be  writ¬ 
ten  as 

An/tfT  -  C/T2  ,  (B-2) 

where  C  Is  the  value  of  the  constant  parameters.  Prom  Equations  (3.1) 
and  (3.2),  assuming  the  time  average  to  be  over  a  short  period  in  which 
temperature  does  not  vary  greatly,  the  time  rate  of  change  of  the  refrac¬ 
tive  index  can  be  given  as 

Dn/DT  -  (An )2/(AT)2  -  cV  .  (B-3) 

Since  the  average  temperature  in  degrees  Kelvin  is  large  with  respect  to 
the  temperature  differences  measured,  (on  the  order  of  tenths  of  degrees 
Kelvin)  a  constant  value  for  Dn/°T  evaluated  at  the  average  air  tempera¬ 
ture  can  be  substituted  in  Equation  (3.3). 

For  the  actual  calculation  of  D  /D_  a  more  precise  formula  for 

D  1 

atmospheric  index  of  refraction  credited  to  Berrell  and  Sears  was  used. 

This  equation  Includes  corrections  for  vapor  pressure  in  the  atmosphere 

and  hither  order  wavelength  factors.  In  order  to  evaluate  the  difference 

quotient  for  the  average  temperature,  the  Barrell  and  Sears  formula  was 

evaluated  for  the  average  temperature  plus  and  minus  one  degree.  The 

average  of  the  two  differences  with  respect  to  average  temperature  was 

then  used  for  D  /Dm  in  the  C  calculations, 
n  T  n 

# 

Smithsonian  Neteorologioal  Tctolss.  6th  sd.,  Washington,  D.C., 

Smithsonian  Institute,  p.  389. 
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The  refractive  index  structure  coefficient,  ru,  is  an  important  parameter  for  optical 
propagation  in  the  atmosphere.  A  theoretical  discussion  of  structure  functions  is 
presented  to  show  the  origin  of  the  refractive  index  and  temperature  structure  func¬ 
tions.  Hie  equipment  used  to  measure  fluctuations  in  the  atmosphere  is  discussed,  and 
the  calculation  of  the  refractive  index  structure  coefficient  from  these  thermal 
measurements  is  shown.  The  details  of  the  experimental  evaluation  of  Cr  and  the 

results  of  C  calculations  are  also  included, 
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